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ABSTRACT 

We present the results of modeling dust spectral energy distributions (SEDs) across the Small 
Magellanic Cloud (SMC) with the aim of mapping the distribution of polycyclic aromatic hydrocarbons 
(PAHs) in a low-metallicity environment. Using Spitzer Survey of the SMC (S'^MC) photometry from 
3.6 to 160 /im over the main star-forming regions of the Wing and Bar of the SMC along with 
spectral mapping observations from 5 to 38 /im from the Spitzer Spectroscopic Survey of the Small 
Magellanic Cloud (S^'MC) in selected regions, we model the dust spectral energy distribution and 
emission spectrum to determine the fraction of dust in PAHs across the SMC. We use the regions 
of overlaping photometry and spectroscopy to test the reliability of the PAH fraction as determined 
from SED fits alone. The PAH fraction in the SMC is low compared to the Milky Way and variable- 
with relatively high fractions (<zpah~ 1 — 2%) in molecular clouds and low fractions in the diffuse 
ISM (average ((/pah) = 0.6%). We use the map of PAH fraction across the SMC to test a number 
of ideas regarding the production, destruction and processing of PAHs in the ISM. We find weak or 
no correlation between the PAH fraction and the distribution of carbon AGB stars, the location of 
supergiant H I shells and young supernova remnants, and the turbulent Mach number. We find that 
the PAH fraction is correlated with CO intensity, peaks in the dust surface density and the molecular 
gas surface density as determined from 160 /xm emission. The PAH fraction is high in regions of active 
star-formation, as predicted by its correlation with molecular gas, but is supressed in H II regions. 
Because the PAH fraction in the diffuse ISM is generally very low-in accordance with previous work 
on modeling the integrated SED of the SMC-and the PAH fraction is relatively high in molecular 
regions, we suggest that PAHs are destroyed in the diffuse ISM of the SMC and/or PAHs are forming 
in molecular clouds. We discuss the implications of these observations for our understanding of the 
PAH life cycle, particularly in low-metallicity and/or primordial galaxies. 

Subject headings: dust, extinction — infrared: ISM — Magellanic Clouds 



1. INTRODUCTION 

Polycyclic Aromatic Hydrocarbons (PAHs) are 
thought to be the carr ier of the ubiquitously o bserved 
mid-IR emission bands (lAllamandola et aI]ll989L among 
others). The bands are the result of vibrational de- 
excitation of the PAH skeleton through bending and 
stretching modes of C-H and C-C bonds after the 
absorption of a UV photon. The emission in these 
bands can be very bright and c an compr ise a significant 
fraction, up to 10-20% (Sm ith et all i2007), of the 
total infrared emission from a galaxy. For this reason, 
PAH emission has been suggested to be a useful tracer 
of the star-formatio n rate, even out to high redshifts 
(jCalzetti et al.ll2007| ). Making use of PAHs as a tracer, 
however, requires understanding how the abundance 
and emission from PAHs depends on galaxy properties 
such as metallicity and star-formation history. 

PAHs also play a number of important roles in the 
interstellar medium (ISM). In particular, these small 
dust grains can domin ate photoelectric heating rates 
(jBakes fc TielensI HOOl ). In dense clouds, PAHs can 
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alter chemical reaction networks by providing a neu- 
tralization route for ionized spe cies (|Bakes fc TielensI 
Il998t iWeingartner fc Draind [200ll ) and contribute large 
amounts of surface area for chemical reactions that occur 
on grain surfaces. PAHs are a crucial component of inter- 
stellar dust so we would like to understand the processes 
that govern their abundance and physical state. 

The life-cycle of PAHs, however, is not yet well un- 
derstood. PAHs are thought to form in the carbon- 
rich atmospheres of s ome evolved stars (jLatteii 119911 : 
ICherchne ff et al.lll992[ ). Emission from PAHs has been 
observed from carb on-rich asymptotic giant branch stars 
oan et all I2007D and more frequently in carbon-rich 
post- AGB stars where the radiatio n field is more e ffec- 
tive at exciting the mid-IR bands (|Buss et al.lll993D . A 
"Stardust" origin (i.e. formation in the atmospheres of 
evolved stars) for the majority of PAH material is contro- 
versial, however, because it has yet to be demonstrated 
that PAHs can be produced in AGB stars faster than 
they are des troyed in the ISM (for a recent review, see 
lDrainel[2009() , i.e. the timescale for destruction of dust 
by SNe shocks is shorter than the timescale over which 
the ISM is enriched with dust from AGB stars (for exam- 
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ple. lJones et al.|[T994D . In addition to destruction by su- 
pernova shocks, PAH material may be destroyed by UV 
fields, a process that can dominate near a hot star or in 
the ISM of low metallicity and/or primordial galaxies. If 
PAHs are mostly not "stardust" , they must have formed 
in the ISM itself, by some mechanism which is not yet 
charact erized. A variety of mechanisms have been sug- 
geste d (iTielens et alj 119871: iPuget fc Legeii fTOSQUHerbsti 
I199H IGreenberg et all 120001) " however there is little ob- 
servational support for any one model as of yet. 

In recent years, observations with ISO and Spitzer 
have allowed us to study the abundance and physi- 
cal state of PAHs in a variety of ISM conditions be- 
yond those we observe in the Milky Way. One of the 
most striking results is the abrupt change in the frac- 
tion of dust in PAHs as a function of metallicity. A 
deficit of PAH emission from low metallicity galaxies has 



been widely ob served (iMa.dderLi.2000 : Engelb racht et al 
2005; Madden et al.ll2006l: iWu et al.ll2 006: Ja ckson et al 



~006i .Engelbracht et al.ll2008D . lEngel'bracht et al.l (120050 



found that the ratio of the 8 to 24 surface bright- 
ness undergoes a transition from a SED with typical 
PAH emission to an SED essentially devoid of PAH 
emission at a metallicity of 12 + log(0/H) ~ 8. The 
weakness of PAH emission in low-me tallicity galaxies 
has been confirmed sp ectroscopically (|Wu et al.l 120061 : 
lEngelbracht et al."2008'). Using the SINGS galaxy sam- 
ple iDraine et a l. (2007) modeled the integrated SEDs 
and determined that the deficit of PAH emission cor- 
responds to a decrease in the PAH fraction rather than 
a change in excitation of the PAHs. They found that 
9PAH (defined as the fraction of the total dust mass that 
is contributed by PAHs containing less than 10^ carbon 
atoms) changes from a median of ^ 4% (com parable to 
the Milky Way PAH fraction of 4.6%: iLi fc P'rainc 2001. ) 
in galaxies with 12 -I- log(0/H) > 8.1 to a median of 
^ 1% in more metal poor galaxies. iMufioz-Mateos et al.l 
(|2009| ) have investigated the radial variation of gpAH and 
metallicit y in the SING S sample and find results consis- 
tent with lDraine et al.l (|2007[ ). 

There have been a number of suggestions as to what 
in the PAH life-cycle change s at low metalhcity le ading 
to the observed deficiency. iGalliano et al.l (|2008| ) sug- 
gested that the delay between enrichment of the ISM 
by supernova-produced dust relative to that from AGB 
stars could lead to a lower PAH fraction at low metallic- 
ity. This model relies on the assumption that supernovae 
contribute a significant amount of dust to the ISM, an 
assumption which is controversial ('Moselev et al .l 119891: 
Drmnc et al. 2003; Krausc et al. 2004; Sugerma n et al.l 
2006t iMeikle et al.l l2007l : iDrainel 12009 ). as weU as long 



timescales for dust production in carbon-rich AGB 
stars, which may be shorter than prev iously thought 
Sloan et al.l 120091) . Fundamentally, the IGalliano et al.l 
20081 ) model assumes a "stardust" origin for PAHs, 



which may not be the case. Other models explaining the 
low metallicity deficiency rely on enhanced destruction 
of PAHs. This can be accomplished through mo re effi- 
cient destruction via supernova shocks (|0'Halloran et al.l 
I2OO6I) or via the harder and more intense UV fields in 
these galaxies (e.g. iMadden et al.l l2006t iGordon et al.l 
[2008h . 

In order to investigate the PAH life cycle at low metal- 



licity, we performed two surveys of the Small Magellanic 
Cloud (SMC) with the Spitzer Space Telescope. The 
SMC is a nearby dwarf irregular galaxy that is cur- 
rently interacting with the MW and th e Large Mag- 
ellanic Cloud. Its proximity (61 kpc; iHilditch et al.l 
, 2005). low metaUicitv (12 -I- log(0/H) - 8, Z ~ 0.2^0; 
iKurt fc DufouHll998l) and tidally disrupted ISM make it 
an ideal location in which to study the life cycle of PAHs 
in an environment very different from the Milky Way. 
The SMC has a low dust-to-gas ratio. ~ 10 times smaller 
thaii in the Milky Way (Bot et aLl 120041 : Eeroy et al] 
|2007| ). leading to more pervasive UV fields. Because of its 
proximity we can observe the ISM at high spatial resolu- 
tion and sensitivity in order to characterize the processes 
driving the PAH fraction. 

Th e PAH fract i on in the SMC has been controver- 
sial. ILi fc Draint^ ([2002^, using IRAS and COBE data, 
found that the PAHs in the SMC Bar contained only 
0.4% of the in terstellar carbon , corresponding to (7pah~ 
0.2%, whereas iBot et al.l (|2004D concluded that PAHs ac- 
counted for 4.8% of the total dust mass in the diffuse 
ISM of the SMC, similar to the Milky Way. In the fol- 
lowing, we present results on the fraction of PAHs in the 
SMC using observations from the Spitzer Survey of the 
Small Magellanic Cloud (S'^MC). We use spectroscopy 
in the regions covered by the Spitzer Spectroscopic Sur- 
vey of the Small Magellanic Cloud (S'^MC) to verify that 
our models for the photometry are correctly gauging the 
PAH fraction. We defer a detailed analysis of the spec- 
troscopy to an upcoming paper (Sandstrom et al. 2010, 
in prep) In Section [5] we describe the observations and 
data reduction, particularly focusing on the foreground 
subtraction and cross-calibration of the IRAC, MIPS and 
IRS observations. In Section [3] we descri be the SED fit - 
ting procedure using the models of Drai ne fc Lil (|2007f ) 
and the modifications necessary to incorporate the S^MC 
spectroscopy into the fit. In Sections [4] and [5] we present 
the results of the SED modeling and discuss their impli- 
cations for our understanding of the PAH life cycle both 
in low metallicity galaxies and in the Milky Way. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spitzer Survey of the Small Magellanic Cloud 
(S^MC) Observations 

We mapped the main star-forming areas of the Bar 
and Wing of the SMC using the IRAC and MIPS in- 
struments as part of the S'^MC project (GO 3316). A 
more comprehensive descriptio n of the observations and 
data reduction can be found in lBolatto et al.l ()2007[ ) and 
iLerov et al.l(|2007l ). The region where the coverage of the 
IRAC and MIPS observations overlap is shown in Fig- 
ure [1] overlayed on the 24 fim image. The mosaics were 
constructed using the MOPEX software provided by the 
SSCQ. The IRAC and MIPS mosai cs were correct e d for^ a 
number of artifacts as described in lBolatto et al.l ()2007D . 
The most important of these for the purposes of this work 
is the large additive gradients at IRAC wavelengths (pri- 
marily 5.8 and 8.0 /im) caused by varying offsets in the 
detectors and the mosaicing algorithm implemented in 
MOPEX. We will discuss these gradients further in Sec- 
tion l2.3l since they become important for determining the 
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TABLE 1 

S'^MC Observation Details 



Band 


Map Noise Level 


CaL Uncertainty 


Resolution 




(MJy sr-i) 


(%) 


(") 


3.6 


0.015 


10=" 


1.66 


4.5 


0.017 


10=" 


1.72 


5.8 


0.055 


10=" 


1.88 


8.0 


0.042 


10=" 


1.98 


24 


0.047 


4 


6.0 


70 


0.664 


7 


18 


160 


0.695 


12 


40 



Note. — The noise levels for the IRAC maps shown here have 
not been multiplied by the extended source calibration. 
^ Calibration uncertainties in the IRAC bands have been in- 
creased to 10% because of the extended source corrections. 



foregrounds present in our maps. 

Since the initial processing of the MIPS mosaics, the 
cahbration factors recommended by the SSC have been 
revised. We correct the mosaics to use the recommended 
factors of 0.0454, 702 and 41.7 MJy sr^^ per instrumen- 
tal data unit, which are 3%, 11% and 0.7% diff erent 
from the earlier va lues used bv lBolatto et al.l ()2007[ ) and 
ILerov et all (|2007D . In addition, the 70 /im observations 
suffer fr om non-linearities at high surface brightness as 
noted bv lDale et al.l (|2007l ). Although relatively few pix- 
els in our map are affected by the correction, these re- 
gions in particular are most likely to overlap with our 
spectroscopic observations. We use the most recent non- 
linearity correction described in Gordon et al. (2010), in 
prep. No calibration correction for extended emission is 
necessary for the MIPS mosaics (Cohen 2009). The la 
sensitivities of the observations are listed in Table [TJ The 
noise level at 70 fim is higher than the predicted detec- 
tor noise due to pattern noise, visible as striping in the 
map (for further dis cussion of 70 /im noise properties, see 
iBolatto et aTll2007| ). Aside from the non-linearity correc- 
tions, we use the MIPS calibr ation factors listed on the 
SSC website, but we note that ILerov et al.l (|2007D found 
an offset between the MIPS 160 fxm and the DIRBE 140 
fiui photometry of the SMC which may be the results of 
a calibration difference. They adjusted the calibration 
by a factor of 1.25 to match DIRBE. We do not apply 
this correction, and we briefly discuss the implications of 
that choice in Section [X^ 

Because the calibration of the IRAC bands is based on 
stellar point sources and we are dealing with extended 
objects, we also apply the extended source calibration 
factors of 0.955, 0.937, 0.772 and 0.73 7 to the 3.6 4.5, 5 8 
and 8.0 bands as recom mended by Rea ch et al.l (j2005D . 
Recent work bv lCoheiTet al. (2007) verified the 36% cor- 
rection at 8.0 /Ltm. The correction factors depend on 
the structure of the emission for each region of the map, 
which ranges from diffuse to point-like. Thus, applying 
on uniform correction factor across the entire map in- 
troduces some systematic uncertainty in our IRAC flux 
densities. To account for these uncertainties we assume 
a 10% calibration uncertainty for the IRAC bands. The 
details of the calibration are discussed further in an Ap- 
pendix and the convolution and alignment to a common 
resolution will be discussed in Section r2. 4.21 

2.2. Spitzer Spectroscopic Survey of the Small 
Magellanic Cloud (S^MC) 



In order to directly probe the physical state of PAHs 
in the SMC, we performed spectral mapping of six star- 
forming regions using the low spectral resolution orders 
of IRS on Spitzer (GO 30491). These observations are 
primarily intended to investigate spectral variations in 
the PAH emission, to be described in Sandstrom et al. 
(2010), in prep. The coverage of the maps are shown 
overlayed on the MIPS 24 /im image from our S'^MC 
observations in Figure [T] 

The spectral coverage of the low-resolution orders of 
IRS extends from 5.2 to 38.0 //m, covering the major 
PAH emission bands in the mid-infrared except the 3.3 
/im feature. The maps are fully sampled by stepping 
perpendicular to the slit one-half slit width between each 
slit position (1.85"and 5.08", for SL and LL respectively). 
The LL maps are made of 98 pointings perpendcular to 
the slit and 7 pointings parallel for a coverage of 493" x 
474", except for the map of N 76 which covers 75 by 
6 pointings and an area of 376" x 395". The SL maps 
are made of 120 pointings perpendcular and 5 pointings 
parallel covering an area of 220" x 208" in each map, 
except for the region around SMC Bl where we use 60 
by 4 pointings covering an area of 109" x 156". 

The spectra were processed with either version 15.3.0 
or 16.1.0 of the IRS Pipeline. The only difference of 
note between these versions is a change in the process- 
ing of radiation hits that does not affect the results we 
pres ent. The details of the data reduction are discussed 
in iSandstrom et al.l ()2009D . In brief, the maps were as- 
sembled using CubisrrQ, wherein a "slit loss correction 
function" was applied, analogous to the extended source 
correction for the IRAC bands, to adjust the calibra- 
tion f rom point sources to extended objects (Smit h et al.l 
[20071 ). Each mapping observation was followed or pre- 
ceded by an observation of a designated "off" position 
at R.A. 1^9^40" and Dec -73°31'30". This position was 
seen to have minimal SMC emission in our MIPS obser- 
vations. The "off" spectra, in addition to subtracting 
the zeroth-order foregrounds, help to mitigate the effects 
of rogue pixels. Additional bad pixel removal was done 
within Cubism. Beyond 35 /im, there are increasing 
numbers of "hot" pixels which degrade the sensitivity of 
our spectra, we trim the LL orders to 35 /<m to avoid 
issues with the long wavelength data. 

Prior to further reduction steps we determine correc- 
tions to match the SL2 (5.2-7.6 /im), SLl (7.5-14.5 
/im), LL2 (14.5-20.75 /im) and LLl (20.5-38.5 /im) 
orders in their overlap regions. For the SL2/SL1 and 
LL2/LL1 overlap, we find that the offsets are best ex- 
plained by small additive effects, which may be due to 
a temporally and spatially varying dark current analo- 
gous to the "dark settle" effect seen in high-resolution 
IRS spectroscopy. Further discussion of these offsets will 
be presented in Sandstrom et al (2010), in prep. In brief, 
we apply correction factors determined by examining the 
unilluminated parts of the IRS detector, very similar to 
what is done in the "darksettle" software for LH available 
through the SSC. After applying these correction factors 
the orders generally match-up to within their respective 
errors. In regions of low surface brightness small resid- 
ual effects have the appearance of a bump around 20 /im 

^ |http: //ssc ■ spitzer ■ caltech. edu/archanaly/contrlbuted/ cubism/ | 
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Fig. 1. — The coverage of the S^MC and S*MC surveys overlayed on the MIPS 24 fj,m map. The color scale is logarithmic, with the 
stretch illustrated in the colorbar. The red boxes show the coverage of the LLl order maps (the LL2 maps are shifted by ~ 3') and the 
green boxes show the coverage of the SLl order maps (the SL2 maps are shifted by ~ 1'). We also identify the various regions of the galaxy 
by the names we will refer to in the remainder of this paper. 



in the stitched spectra. These offsets only affect a small 
portion of the spectrum and do not measurably alter the 
results of the fit. 

2.3. Foreground Subtraction 
2.3.1. IRAC and MIPS Foreground Subtraction 

In the mid- and far-infrared there are three major fore- 
ground/background contributions that contaminate our 
observations of the SMC: the zodiacal emission, Milky 
Way cirrus emission and the cosmic infrared background 
(CIB). In addition, for the IRAC mosaics there is a pla- 
nar offset introduced by varying det ector offsets and the 
mosaicing algorithm in MOPEX (see iBolatto et~alll2007l 
for more details). These need to be removed from the 
maps to isolate the emission from the SMC. In the fol- 
lowing wc will briefly discuss the foreground subtraction 
that has been carried out on our data and the major 
uncertainties in this process. The approach we use to 
subtract these foregrounds is motivated by the following 
limitations of our observations: (1) at longer wavelengths 
the mosaics do not extend to areas with no SMC emis- 
sion, (2) the observations of neutral hydrogen which we 
use to subtract the Milky Way cirrus emission have a 
lower angular resolution than our Spitzer maps and (3) 
the residual mosaicing gradients in the 5.8, 8.0 and to a 
lesser degree 4.5 fim IRAC bands interfere with directly 
fitting and subtracting the zodiacal light. 



Our foreground subtraction has three steps: (1) we de- 
termine the coefficients of a planar surface that describe 
the combination of zodiacal light and the residual mo- 
saicing offsets, (2) we subtract this plane from each map 
at its native resolution, (3) we convolve the maps to the 
MIPS 160 fim resolution of 40" and subtract the MW 
cirrus foreground and CIB. Throughout this procedure 
we fix the CIB level at 70 and 160 fj,m to be 0.23 and 
1.28 MJy sr~^ as determined from the Spitzer Obser- 
vation Planning Tool (SP0Tj3- We also assume a fixed 
proportionality between the infrared cirrus emission and 
the column density of MW H I (Boulangcr ct al. 1996), 
using coefficients derived from the model of lDraine fc O 
(|2007f ) for MW d ust heated by the lo cal interstellar ra- 
diation field. The lDraine fc Lil ()2007f) mod el reproduces 
the DIRBE observations of the cirrus from lArendt et ahl 
(1998) to within their quoted 20% uncertainties. We con- 
volve the Draine & Li (2007) model emissivity spectrum 
with the IRAC and MIPS spectral response curves to ob- 
tain the coefficients, which are listed in Table [2] We use 
a map of Galactic hydrogen ob tained by combining data 
from ATCA and Parkes (see iStanimirovic et aLl 119991 
for further information on the technique and the original 
observations) re-reduced and provided to us by E. MuUer 
(private communication) to subtract the MW cirrus con- 

^ |http: //ssc . spitzer . caltech.edu/docunients/background/bgdoc_releaseTi 
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TABLE 2 

Milky Way Foreground 
Coefficients 



Band 


Coefficient 




(MJy sr-i (lO^l H)"!) 


3.6 


0.018 


4.5 


0.006 


5.8 


0.070 


8.0 


0.215 


24 


0.162 


70 


2.669 


160 


11.097 



tribution. 

The first step in our foreground subtraction is to deter- 
mine the coefficients describing the planar contributions 
from the zodiacal light and the mosaicing offsets. Over 
the area of the SMC, the zodiacal light is well-described 
by a plane and the gradients introduced by the mosaicing 
algor ithm also have a planar dependence (jBolatto et alJ 
l2007f ). Unfortunately, there are no "off" locations where 
we could simply fit a plane to the foreground level, since 
no regions are totally free of SMC emission in the S'^MC 
maps. At each position in the map the surface brightness 
is a combination of: SMC emission, MW cirrus emis- 
sion, zodiacal light, mosaicing offsets and CIB. To iso- 
late the planar component, we first extract photometry 
for a number of 200" x 200" regions from the IRAC and 
MIPS maps. We choose the regions to be outside of star- 
forming areas, where the SMC emission is dominantly 
from dust heated by the general interstellar radiation 
field, where we can assume a proportionality between the 
dust emission and the SMC H I column. We choose the 
region size of 200" x 200" to be larger than the 98" beam 
of the H I observations and large enough to robustly 
determine the mean level in each box even when con- 
taminated by point sources at the shorter wavelengths. 
We then subtract off the Milky Way cirrus contribution 
for those regions determined with the coefficients in Ta- 
ble [1 and the MW H I map, and, for the 70 and 160 
/im maps, the CIB level. We call this MW cirrus and 
CIB subtracted value S,y_rcsid- Sj/,rosid is a combination of 
the zodiacal light, emission from the diffuse ISM of the 
SMC and whatever residual gradients there remain from 
the mosaicing for the IRA C maps. Next, usin g the S MC 
neutral H I observations of lStanimirovic et al.l ()1999[ ) , we 
perform a least-squares fit to the foreground values with 
the following function: 

= A{l + BAa + CAs) + D X HIsMC- (1) 

Here A, B and C are the coefficients describing a plane; 
Aq and are the gradients in Right Ascension and 
Declination across the SMC; D is dust emission per H in 
the given waveband; and Si/^resid is the residual emission 
in each box after subtracting the MW foreground and 
CIB components. With this fit we determine the coeffi- 
cients of the best fit planar foreground while excluding 
emission in the map coming from the SMC itself. 

We take two additional steps to improve the determi- 
nation of the planar foreground coefficients: (1) we use 
the 2MASS point source catalog to avoid regions of high 
stellar density in the IRAC bands so we do not subtract 
unresolved starlight, which can masquerade as a fore- 



ground, and (2) we fix the gradient of the zodiacal light 
[B and C) at 70 and 160 /zm to the results of the fit at 24 
/zm. Since there are no mosiacing offsets for the MIPS 
bands, the zodiacal light is by far the dominant fore- 
ground at 24 ^Tci and the zodiacal light should have the 
same spatial dependence at all of the MIPS wavelengths, 
fixing the zodiacal light gradient to what we measure at 
24 /j,m is more effective than trying to fit for those coeffi- 
cients at the longer wavelengths. The fixed pattern noise 
at 70 /im and the increasingly dominant SMC emission 
at 70 and 160 /im make the fitting procedure less robust 
compared to the 24 /j,m results. 

The results of the fits are listed in Table [31 For the 
IRAC bands at 3.6 and 4.5 /j-m we do not detect any 
emission correlated with the SMC neutral hydrogen. At 
4.5 /im, there is a quite significant residual gradient from 
the mosaicing procedure. At 3.6 /^m the zodiacal fore- 
ground and its gradient are not detected at the sensitivity 
of the map. At both 3.6 and 4.5 unresolved starlight can 
play a role in the foreground determination, so we choose 
regions to avoid high stellar densities. In Table [3] we also 
list the standard deviations of the fit residuals to show 
the quality of the foreground determination for the mo- 
saics. Finally, we subtract the planar fit listed in Table[3] 
from the mosaics at their full resolution. 

Because it is at lower resolution (98"), we subtract the 
Milky Way foreground after convolving to the MIPS 160 
/im resolution 40"). At the signal-to- noise ratio of the 
map, the Galactic H I in front of the SMC does not show 
any high contrast features, although there is a gradient 
ranging between 2 and 4x10^° cm~^ across the region. 
We note that the resolution difference between our mo- 
saics and the H I observations means that there could be 
features with spatial scales less than 98" that are not ad- 
equately subtracted from our maps. However, assuming 
a distance of 1 kpc for the H I foreground of the Galaxy, 
features left in our map would have to have spatial scales 
of ^ 0.5 pc. The typical foreground of Milky Way gas 
in this direction is ~ 3 x 10^*^ cm~^, an 0.5 pc cloud of 
cold neutral medium having a typic al volume density of 
- 40 cm-3 (|Heiles fc Troland I [2001 ) would contribute a 
column of 6 X 10^^ cm~^, a factor of 5 less than the aver- 
age foreground we see. Thus, we expect that inadequate 
subtraction of small scale structure will not greatly affect 
the quality of the foreground subtraction. 

2.3.2. IRS Foreground Subtraction 

The "off" observations for each spectral cube will re- 
move foreground emission to first order. However, there 
are gradients between the "off" position and the cubes 
that are still present in the data. To remove the remain- 
ing foregrounds we add back to the spectra the difference 
between the "off" position and the map position at each 
wavelength calculated using the zodiacal light predictions 
from SPOT a nd th e MW cirrus emission spectrum from 
iDraine fc L? (2007) multiplied by the MW H I column 
density at those locations. 

2.4. Further Processing 

2.4.1. Point Source Removal 

The SED models which we use to determine the PAH 
fraction assume a stellar component in the Rayleigh- 
Jeans tail with (xv^. In cases where this condition is 
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TABLE 3 

Foreground Properties 



Band 


A 


B 


C 


D 


Std. Dev. of Fit 




(MJy sr-i) 


(dei 




(MJy sr-i (10^0 H)"!) 


(MJy sr-i) 


3.6 


<0.004 








0.011 


4.5 


-0.11 


0.019 


4.074=^ 




0.014 


5.8 


4.99 


-0.002 


0.002 


9.14 X lO-"* 


0.023 


8.0 


3.69 


0.0014 


0.0014 


7.43 X 10-5 


0.023 


24 


22.34 


-0.0026 


0.0010 


1.63 X 10"^ 


0.039 


70 


6.43 


-0.0026'' 


O.OOIO"" 


3.97 X 10-2 


1.23 


160 


1.83 


-0.0026*= 


0.0010*= 


1.81 X 10-^ 


3.30 


Note. 


— Sec Scctionl2.3. 11 for a 


description 


of the coefficients. 





^ Large gradient at 4.5 /^m due to detector offsets and mosaicing algorithm. 
^ Fixed from fit at 24 /xm. 



not met, the presence of a stellar point source can corrupt 
the results of our fitting. In the vast majority of cases we 
see no problems relating to stellar point sources in the 
SMC, and thus we do not perform a comprehensive point 
source extraction on the observations. In addition, we are 
not able to remove any unresolved stellar contribution, so 
obtaining a map with the stellar component completely 
removed would require more detailed modeling. Instead, 
we mask out bright point sources which do not have a 
well behaved SED in the mid-IR. These sources typically 
fall into one of two categories: very bright stars which 
are saturated at one or more of the IRAC bands or stars 
with non-typical infrared SEDs, such as YSOs or carbon 
stars. We mask these objects out with a circular aperture 
which we fill in with the local background value. 

2.4.2. Convolution and Alignment 

The IRAC and MIPS mosaics are all convolved to 
match the resolution o f the 160 /xm o bserva tions using 
the kernels derived by iGordon et alJ (|2008l ) and avail- 
able from the SSC. After convolution we regrid the maps 
to match the astrometry and pixel scale of the 160 fim 
mosaic. 

For the IRS cubes, the convolution and alignment pro- 
cedure involves a few additional steps. We convolve di- 
rectly with the 160 ^m PSF at each wavelength, which is 
much larger than the PSF of IRS even at its long wave- 
length end, so wavelength dependence of the PSF makes 
little difference to the final map. After convolution we 
align the cube s using the p olygon clipping technique de- 
scribed in Sand strom et al.l (;2009) . Throughout this pro- 
cess, we appropriately propagate the uncertainty cubes 
produced by Cubism. Although the AORs are the largest 
possible size given the observation length limitations on 
Spitzer, the maps are only at most a few resolution el- 
ements across at 160 /im. Thus, we take some care to 
make sure that the spectra we use are only those where 
the PSF is not sampling regions outside the observed 
cube. For the following analysis we use only spectra 
where 90% of the area of the PSF or more is within 
the observed cube, which cuts down the number of vi- 
able spectra we can extract from the cubes to 63. The 
small number of resolution elements across the cube re- 
sults in the introduction of scatter into the comparison 
between the MIPS and IRAC mosaics and the spectral 
maps. This is the result of information from outside the 
cube not being "convolved in" . This scatter represents a 
fundamental limitation of our observations. We estimate 



the magnitude of this scatter to be ~ 10% by comparing 
convolved and aligned IRAC 8.0 and MIPS 24 /im maps 
that have been cropped to the coverage of the spectral 
cubes to those that have not. For the SMC Bl cube we 
do not convolve or align the map because of the small 
number of resolution elements and the loss of signal due 
to the subsequent regridding to match at 160 /im. In- 
stead we extract the spectrum for the whole SMC Bl 
spectral map and extract matching photometry from the 
IRAC and MIPS observations. After convolution and 
alignment we apply the previously mentioned correction 
factors to match the orders and stitch the cubes together 
in the overlaping spectral regions. 

3. SED MODEL FITTING 

We use the SED models of lDraine fcLl (|2007t ) to deter- 
mine the dust mass, radiation field properties and PAH 
fraction using our MIPS and IRAC mosaics in every in- 
dependent pixel of the map where all of the MIPS and 
IRAC measurements are detected above 3a. We also per- 
form simultaneous fits to photometry and spectroscopy 
in the regions covered by S^MC as described in Sec- 
tion l3.3l To distinguish between these two types of fits we 
introduce the following terminology: "photofit" refers to 
the best-fit model using only the photometry and "pho- 
tospectrofit" refers to the best-fit model to the combined 
photometry and spectroscopy. 

The model fit involves searching through a prc-made 
grid of models and finding the model which minimizes 
the following pseudo-x^: 



X 



E 



(-Fobs.fc — {Fu , mo del 



(2) 



^obs,& 



model, 6 



Here Fobs,& is the observed surface brightness in band 
b, (^!y,modoi)6 IS the modcl spectrum convolved with the 
spectral response curve of band b, Uobs.b is the uncer- 
tainty in the observed surface brightness and a^ode\,b is 
a factor which allows us to account for th e systematics 
asso ciated with the modeling. Following iDraine et al.l 

(|2007;) we use 0-model,6 = 0.1i^,niodel,fc- 

The radiation field heating the dust is described by a 
power-law plus a delta function at the lowest radiation 
field: 



dMc 
dU 



= (l-7)AfD(5(f/-C/™„)+7Afi 
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TABLE 4 



IDraine & 


Lil 1120071) 


Model Parameters 


Parameter 


Range 




> 




> 


gpAH 


0.4-4.6% 


^ min 


0.6-30 


Umax 


10^ - lO'^ 


7 


> 


o 


1.5-2.5 



Here U is the radiation field in units of the i ntersteUar ra- 
diatio n field in the Solar neighborhood from lMathis et ahl 
(fT98l . 7 is the fraction of the dust heated by the min- 
imum radiation field, a is the power-law index of the 
radiation field distribution and is the dust mass sur- 
face density. This parametrization allows us to approxi- 
mately account for both the dust heated by the general 
interstellar radiation field and the dust heated by nearby 
massive stars. 

The adjustable parameters in the model are: the stellar 
luminosity per unit area (i^*), the dust mass surface den- 
sity (M^)), the PAH fraction ((Jpah), the minimum and 
maximum radiation field (Umin and Vmax), the fraction 
of dust heated by the minimum radiation field (1 — 7) 
and the exponent of the radiation field power- law {a). 
r2*, M/) and 7 are continuous variables which normalize 
the grid of models to match the observations. The other 
variables are adjusted in discrete steps and can have the 
range of values listed in Table S] Using the results of the 
fit we compute a number of useful parameters describ- 
ing the outcome. These are fpoR, the fraction of the 
total infrared power radiated by dust grains illuminated 
by radiation fields U > 10^, and U, the average radiation 
field. 

In t he analysis we present here, we use the IDraine fc Lil 
(|2007f) Milky Way (MW) dust model, so we take a mo- 
ment to justify and explain this decision and its impli- 
cations. In order to model the dust emission spectrum, 
it is necessary to choose a physical model for the dust 
which prescribes the abundances of different grains and 
their size distributions. At present, a dust grain size dis- 
tribution model specific to the SMC including a variable 
PAH fraction and conforming to constraints on the dust 
mass and raw materials available has not been estab- 
lished, and producing such a model is complex and be- 
yond the scop e of this paper (it would requ ire extending 
the model by iWeingartner fc Draind 1200 ID . Note that 
because we are primarily interested in the PAH fraction, 
we must employ a model that includes a size distribution 
that extends into the PAH regime. In addition, we aim 
to compare our results with other studies of the PAH 
fraction (particularly that on the SINGS s ample) which 
use th e Mil ky Way dust model developed by IDraine fc Lil 
(|2007i ) and IDraine eTall (|2007l ). 

There is evidence that the dust grain size distri- 
butions in the SMC and the Milky Way are differ- 
ent dR odrigucs et ^1 119971: IWeingartner fc Draind [20011 : 
IGordon et al.i |200^ and that low metallicity galaxies 
in general may have an excess of small dust grains rel- 
ative to the Milky Way. These differences, however, 
should not affect our measurements of the PAH fraction. 



IDraine et al.l (|2007D found that the dust mass and PAH 
fraction from the best fit models for the galaxies in the 
SINGS sample that fell in the range of (/pah covered by 
the LMC/SMC dust models were not significantly differ- 
ent if the Milky Way model was used instead, hence they 
employed Milky Way models. Given our desire to find 
the range of PAH fractions in the SMC and the benefit 
of being able to comp are our results directly with those 
of IDraine etld] (|2007f ). we will fit MW dust models to 
the photometry. 

We expect the choice of the MW dust models will intro- 
duce some systematic effects into the results of the model 
fitting. In particular, there is evidence that the SMC, like 
other low- metallicity galaxies, has a larger contribution 
from small grains compared to the Milky Way which re- 
sults in "excess" emission at 24 and 70 iim ( Gordon ct afl 
200 llBot et"an[200l iGalhano et aDIFoOSt iBernard et aLl 
20081). The effect of such a population on our modeling 
will be to increase the best-fit radiation field in order 
to match the 24 and 70 fj,m brightness. Since we use a 
power-law distribution of radiation fields, and the equi- 
librium emission from large grains determines the min- 
imum radiation field, the general effect is to decrease 
the power-law exponent and increase Umax, making a 
small fraction of the dust be heated by a more intense 
radiation field. Because the PAH emission is produced 
by single-photon heating, gpAH is simply proportional to 
the ratio of the PAH emission in the 8 /im band to the 
total far-IR emission, which is robustly constrained by 
the 70 and 160 ^m photometry. Hence, the derived gpAH 
is relatively insensitive to variations in the Umm, ^max 
and a, provided the observed total infrared luminosity is 
reproduced by the models. The regions with overlapping 
spectroscopy will provide a good test for this reasoning, 
since the 5 — 38 /im continuum in conjunction with the 24, 
70 and 160 //m SED provides more stringent constraints 
for the models. We will revisit this subject in Section 

3.1. The IRAC 4-5 l^-m Bra Contribution 

The models we employ only calculate contributions to 
the IRAC and MIPS photometry from starlight, dust 
continuum and dust emission features. However, the 
emission spectrum from the ISM of the SMC is likely 
to contain a number of emission lines, particularly from 
H II regions, that contaminate the photometry. Most 
of these lines are very weak with respect to the dust 
continuum except in the case of the Brackett a hydro- 
gen recombination line at 4.05 /xm. Recent work by 
iSmith fc HancocH (|2009l ) has shown that low-metallicity 
dwarf galaxies with recent star formation can have a Bra 
contribution that is a significant fraction of their inte- 
grated 4.5 /im flux. This problem is exacerbated within 
the H II regions themselves, where the 4.5 /xm emission 
can be alm ost entirely from B ra (see for instance work 
on M 17 bv lPovichet alll2007l) . 

The SED models contain starlight continuum with a 
fixed 3.6 nm to 4.5 fj,m ratio. If there is excess 4.5 
/im emission due to the hydrogen line, the starlight 
continuum will be too high, altering the ratio of stel- 
lar to non-stellar emission at 5.8 and 8.0 /im. To cor- 
rect for the Bra emission, we use the Magellanic Cloud 
Emission Line Survey imag e of Ha (kindly provided by 
C. Sm ith and F. Winkler; 'Smi th fc The MCELS Teainl 
119991) along with the Case B factors at 10,000 K from 
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lOsterbrock fc Ferlandl (|2OO60 to convert Ha to Bra. This 
correction will underestimate the Bra emission where 
there is extinction, but since the intrinsic extinction in 
the SMC is low to begin with we expect this estimate 
of Bra to be adequate. This Bra correction is highest 
in H II regions, where we find Bra emission can account 
for 20-30% of the total 4.5 /im emission. Outside of H II 
regions, the correction is negligible. 

3.2. Additional Systematic Uncertainties 

iLerov et"al] (|2007D investigated the integrated far-IR 
SED of the SMC, using observations from IRAS, ISO, 
DIRBE and TOPHAT. They found that the MIPS 160 
/im photometry is high by approximately 25% compared 
to the values predicted by interpolating from the DIRBE 
observations at 140 and 240 /im, a significant offset given 
how close in wavelength the DIRBE 140 and MIPS 160 
bands are. More recent MIPS observations of the SMC 
obtained through the SAGE-SMC legacy survey (Gordon 
et al 2010, in prep.) agree very well with the S'^MC pho- 
tometry, suggesting that this offset is real. The source 
of this offset is not known, it may be partly due to a 
calibration difference between the instruments or to [C 
II] 158 /im emission contributing in the MIPS bandpass. 
In their analysis, ILerov et all (j2007D divided the S^MC 
map by a factor of 1.25. We do not use this factor in our 
analysis. Dividing the 160 /im photometry by 1.25 would 
produce a lower PAH fraction from our analysis because 
the 70/160 /im ratio would increase, the radiation field 
needed for the large grains to acheive the necessary tem- 
perature would be higher and, given the emission at 8 
/im, the PAH fraction necessary will be lower. 

If the contribution of the [C II] line is significant, there 
may be a systematic o ffset in the PAH fraction we de- 
termine. iRubin et al] ()2009f ) find that [C II] and 8 /im 
emission are correlated in the LMC. If regions of the 
SMC with 8 /im emission have significant [C II] emission 
in the 160 /im band, we will systematically overestimate 
the PAH fraction, since it will seem that the large grains 
are colder than they really are. Future Herschel observa- 
tions of the [C II] line in the SMC will help to understand 
what the level of contamination of the 160 /im observa- 
tion. 

3.3. Simultaneous SED and Spectral Fitting 

In the following, we use the S^MC data is to investigate 
the reliability of the determination of qpAH from SEDs 
alone (i.e. photofit). The photospectrofit model fitting 
procedure is very similar to the SED fits, and involves 
searching the same grid of models for the best fit. 

In using SED fits to measure the PAH fraction, we 
must assume that the IRAC 8.0 /im band, which sam- 
ples the 7.7 /im PAH feature, traces the total PAH emis- 
sion. This has been seen to be a good assui nption when 
looki ng at the integrated spectra of galaxies ([Smith et al.l 
|2007[ ). but variations in the relative band strengths can 
be averaged out on galaxy-scales and the 7.7 /im fea- 
ture may not trace the total PAH emission as effectively 
within individual star-forming regions. The simultane- 
ous SED and spectral coverage provided by S'^MC and 
S^MC allow us to test the effectiveness of SED fits to 
determine (?pah, since the photofit results will be solely 
determined by the 7.7 feature, while the photospectrofit 



results will, to first-order, match the total PAH emission 
in the mid-IR. 

The photospectrofit models will only reproduce to to- 
tal PAH emission to first order because the models use 
fixed spectral profiles for the various PAH bands, and a 
fixed PAH ionization fraction. There are variations in 
the PAH bands observed in the SMC which are not re- 
produced in detail by the models. However, fitting the 
full mid-IR spectrum in the S^MC regions will approx- 
imately reproduce the total PAH emission in all of the 
major mid-IR bands, and therefore be a better tracer 
than the 7.7 feature alone. Although the comparison be- 
tween the photofit and photospectrofit results is not the 
ideal way to quantify the dependence of gpAH on the band 
ratios, it is the best that can currently be done without 
further modeling which is outside the scope of this paper. 

Prior to fitting the IRS spectroscopy, we remove the 
emission lines from the spectrum using th e PAHFIT 
spectral fitting package ([Smith et al.l I2007D . We also 
use the PAHFIT results to estimate the uncertainties 
in the spectra employing the following procedure. The 
IRS pipeline produces uncertainties based on the slope 
fits to individual ramps and propagates these through 
the various reduction steps. These errors do not include 
any systematic effects, and underestimate the scatter in 
our spectra significantly. In addition, the assumption 
that the uncertainties are random, uniform and uncor- 
related for propagation through our analysis does not 
apply: striping in the spectral cube does not average 
out spatially or spectrally. To get a better idea of the 
uncertainties, we estimate the average deviation of the 
spectrum from the best PAHFIT result in a 5 pixel-wide 
sliding window. Although this technique may artificially 
increase the uncertainties in regions that are poorly fit 
by PAHFIT, we find that uncertainties estimated in this 
manner are much more reasonable than those propagated 
from the pipeline. Finally, because the number of points 
associated with the spectrum is much larger than the 
SEDs it was necessary to artificially increase the weight 
of the SED points to achieve a decent fit. We find that 
applying weighting factors of 40 to all the photometric 
points yield reasonable joint fits. 

4. RESULTS 
4.1. Photofit and Photospectrofit qpAn Consistency 

Figure [2] shows the values of gpAH for the best fit mod- 
els in the regions with overlapping photometry and spec- 
troscopy. Because of the gridded nature of the models, 
there are typically a number of points overlapping in the 
plot, so we additionally show a histogram of the values 
on either axis. The majority of points in our spectral 
map regions have qpAK at the lower limit of the range 
and the majority of those points yield the same value of 
qpAH from the photofit and photospectrofit results. At 
the high end of the range of (?pah , we also see good con- 
sistency between the photofit and photospectrofit values. 
In the intermediate regions, the photofit models tend to 
underestimate qpAH by a small amount. Excluding all of 
the points having the minimum PAH fraction, the pho- 
tospectrofit gpAH value is, on average, larger than the 
photofit value by only AgpAH~ 0.23%. 

Despite the good agreement between the photofit and 
photospectrofit gpAH, the best fit models in the two cases 
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show striking systematic differences. Figure 12] shows a 
comparison of a few of the photofit and photospectrofit 
models, chosen to highhght the range of gpAH we mea- 
sure from the spectra. Table [5] shows the best fit radia- 
tion field parameters for the plotted SEDs and spectra. 
In all cases, the spectroscopic information shows that 
the mid-IR continuum below the PAH features is lower 
than the photofit model predicts. The differences arise 
because the photofit models have the observed SED un- 
constrained over the factor of ~ 3 gaps in wavelength be- 
tween 8 and 24 fun and between 24 and 70 /im, whereas 
the photospectrofit models add continuous constraints on 
the SED from 5 to 38 /im. The photofit models tend to 
overpredict the continuum between 8 and 24 ^m, and to 
underpredict the continuum between 24 and 70 /zm, using 
a values that are too small, and leading to overestimation 
of fpDR- The overprediction of the 8 — 24 /im continuum 
leads to underprediction of (?pah, but as seen in Figure [21 
the bias is not large, amounting to AqpAH~ 0.23%, al- 
though in some cases the errors are larger. In general 
the dust surface density and stellar luminosity are not 
changed in a systematic way. What this amounts to is a 
redistribution of the radiation field to increase the frac- 
tion of dust that is being heated by very high radiation 
fields and to decrease the radiation field necessary in the 
diffuse ISM. In fact, the average radiation field in the 
photofit and photospectrofit models is similar. We show 
a series of plots illustrating these changes in Figures |4] 
and[5l 

As previously discussed, the value of qpah is essen- 
tially proportional to the ratio of the power in the 8 /im 
PAH emission features to the total far-IR power, and is 
therefore relatively insensitive to variations in the other 
fitting parameters. The agreement between the two best 
fit gpAH values is a strong indication that the technique is 
robust even though our dust model is not tailored for the 
SMC and does not have variable band ratios. We note 
that the spectroscopic maps are preferentially located in 
star-forming regions, and most cover H II regions. In 
these spots in particular, the radiation field will deviate 
the most from the general interstellar field. Over the rest 
of the SMC, the increase in fpoR will not be as dramatic. 
We also note that the largest differences in qpah from the 
spectoscopic and photometric comparison occur in N 22, 
which contains a point source that is highly saturated at 
24 /im. We have attempted to exclude the regions of N 
22 affected by the saturation, but if there were excess 24 
/im emission from the PSF wings of this source contam- 
inating the photometry, it may artificially drive up the 
PDR fraction and change qpah more drastically. 

Our conclusions from the comparison of the photofit 
and photospectrofit models are that the qpah values are 
in agreement. The radiation field parameters from the 
photospectrofit models reflect a redistribution of the ra- 
diation such that the average stays the same while the 
fraction of dust heated by more intense fields increases 
and the minimum field decreases. These shifts are most 
likely not as large in most regions across the SMC as 
they are in the star-forming regions we probe with spec- 
troscopy. Finally, for regions with intermediate values of 
9pah we recognize the fact that we may underestimate 
the PAH fraction by a few tenths of a percent on av- 
erage, however this small difference does not affect our 
conclusions. 



4.2. Results of the SED Models Across the SMC 

In Figure [6] we show representative SEDs and photofit 
models from four locations in the SMC and in Figure [7] 
we show the gpAH from the photofit models at every pixel 
in our map. All pixels in our map have gpAH less than the 
average Milky Way value (qpah.mw = 4.6%). One of the 
noticable features of this map is the large spatial varia- 
tions in the PAH fraction, from essentially no PAHs to 
approximately half the Milky Way PAH fraction in some 
of the star-forming regions — a range that spans nearly 
an order-of-magnitude. 

The average PAH fraction in the region we mapped is 
(9pah) — 0.6%, determined by the following average: 



(^pah) 



(4) 



Given the minimum value of qpah allowed by our models, 
this value is in good agreement wit h the SW Bar average 
determined bv iLi fc Draind (12002 ) but is 8 times lower 
than the average from iBot et al.l ( 20041 ) . 

Previous studies of low metallicity dwarf galaxies have 
shown large variatio i is in 8 /zm surface b r ightness (e.g. 
Cannon et all l2006t lJax;kson et all l2006t iHunter et ahl 
2006; Wal ter et al.ll2007D which we also identify in the 



SMC. Some of the regions that are brightest at 8 /x m have 
relatively low PAH fractions (c.f. ICannon etall [20061. 
To illustrate, we overlay the (Jpah map on the IRAC 8 
/im mosaic in Figure [S| There are a number of regions 
where 8 /im emission is very bright while the PAH frac- 
tion is low. In particular, N 66 and the Northern region 
of the SW Bar stand out as very bright 8 /im sources 
which have relatively low qpah- A representative spec- 
trum of N 66 is shown in the bottom panel of Figure [31 
illustrating the low qpah in this region. 

To evaluate the use of the 8/24 ratio as an indicator 
of PAH fraction, we show in Figure [HI a plot of qpah vs 
8/24, with the average for each value of qpah overlayed. 
The PAH fraction is correlated with the 8/24 ratio, even 
on small scales in the SMC. However, it is evident that 
the correlation is weak, with a large range in qpAH for a 
given value of the 8/24 ratio. 

The 8/24 r atio spans the region where 
lEngelbracht"etaI[ (2005) see a transition from galaxies 
with evidence for PAH emission to galaxies which show 
no PAH features, which is what one would expect given 
the metalhcity of the SMC. The resuhs of our SED 
modeling indicate that the PAH fraction, though always 
lower than (7pah,mWj is not uniform across the SMC. 
There are regions with qpah within a factor ~ 2 of 
qpA H.MW and regions w here <7pah is at the lower limit of 
the lDraine fc Lil (j2007[) models (~ (7pah,mw/10)- Since 
there are not comparable resolved maps of gpAH in a 
sample of galaxies spanning this transition zone, we 
cannot explain the trend in PAH fractions by looking 
at the SMC alone. However, the global PAH fraction 
we measure for the SMC (((/pah) ^ 0.6%) is driven by 
the very low qpak over the majority of the galaxy and 
the large regional variations in (/pah make it unlikely 
that we are observing a uniform decrease in the SMC 
PAH fraction. // the SMC is typical of galaxies at this 
metallicity, the transition represents a decrease in the 
filling factor of the PAH-rich regions, rather than a 
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Fig. 2. — gpAH from photofit and photospcctrofit models. Because of the gridding of the models, a number of points can overlap on 
this plot. The gray line shows the average offset between the two values excluding the points which have qpAii= 0.4% while the black line 
shows a one-to-one relationship. Most of the points have a best fit value for gpAH at the lower limit of the model range. 



TABLE 5 

Parameters of Selected Photofit and Photospectrofit Models Shown in Figure [3] 



Panel 




Photofit 








Photospectrofit 




a 


Uniin 


Umax 


7 




a 


Umin 


Umax 


7 


a 


1.70 ± 0.08 


3.0 ± 1.0 


7.0 


4.9 ±3.8 X 10' 


4 


1.50 ± 0.70 


2.0 ± 0.8 


3.0 


0.1 ±4.1 X 10^1 


b 


1.50 ± 0.04 


1.5 ± 0.5 


7.0 


2.7 ±2.2 X 10" 


-5 


1.70 ± 0.41 


1.2 ± 0.2 


3.0 


0.2 ±3.4 X 10-1 


c 


1.80 ± 0.06 


5.0 ± 1.1 


7.0 


2.1 ±0.9 X 10" 


-3 


2.30 ± 0.01 


2.0 ± 0.8 


4.0 


4.1 ±2.9 X 10-1 


d 


2.30 ± 0.19 


5.0 ± 2.4 


7.0 


9.0 ±3.0 X 10- 


-1 


2.20 ± 0.29 


3.0 ± 12. 


5.0 


1.0 ±0.7 X 10° 



uniformly low global PAH fraction. 

The SED fits provide a number of parameters describ- 
ing the radiation field. In Figure [101 we show two panels 
which illustrate the average radiation field U and the 
PDR fraction fpoR- Regions where the PDR fraction 
is high tend to correspond to H II regions, as expected. 
We overlay representative contours of Ha on the map of 
fpDR and U to highlight the brightest H II regions in 
the Cloud. We also measure the total dust mass and the 
stellar luminosity at each pixel. We find that the dust 
mass from our fits agrees well with previous results from 
iLerov et all (|2007[) using the same MIPS observations 
despite different methodologies. Figure [TT] s hows a com- 
pariso n of our dust mass results with those of lLerov et al.l 
pool and we find that our masses are lower by ~ 30%, 
we ll within the ^ 50% systematic uncertainties claimed 
bv ILerov et all (|200l . 

4.3. Spatial Variations in the PAH Fraction 

Since we see clear spatial variations in the PAH frac- 
tion, we discuss in the following section what sort of 



conditions foster high PAH fractions in the SMC. We 
observe three trends: (1) the PAH fraction is high in re- 
gions with high dust surface densities and/or molecular 
gas as traced by CO, (2) the PAH fraction is low in the 
diffuse interstellar medium and (3) the PAH fraction is 
depressed in H II regions. 

In Figure [12] we show the qpah map overlayed with 
contours of ^ ^CO (J= 1 — 0) em ission from the NAN- 
TEN survey izuno "ei~all[200l . There is a strong cor- 
relation between the presence of molecular gas and the 
regions with higher PAH fraction. In Figure [T3| we show 
a histogram of gpAH for lincs-of-sight with detected CO 
emission in the NANTEN map of the SMC and those 
without. The mean value of qpah for a line of sight with 
CO is twice that for a line of sight without CO. In ad- 
dition, there are no lines of sight through only atomic 
gas which have gpAH higher than ~ 1%. We note that 
the CO map has much lower resolution than our map of 
9PAH, so the association of PAHs with CO is most likely 
stronger than the evidence we present here. 

The association of PAH emission with star-forming re- 
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Fig. 3. — Four examples of the photofit and photospectrofit models for regions with overlapping photometry and spectroscopy. The 
IRS spectrum is shown in gray and the MIPS and IRAC photometry are shown with diamond-shaped symbols. For clarity we have not 
overplotted the model photometry. The photofit model is shown as a dashed black line and the photospectrofit model is shown as a solid 
black line. The differences between the models represent a redistribution of the radiation field, the parameters of which are listed in TablelSl 
The R.A. and Dec position of each of these regions is listed in the upper left corner of the plot. Panel d) shows a representative spectrum 
from the N 66 region. 



gions and molecular clouds versus the diffuse ISM of a 
galaxy is a matter of debate, and may vary depending on 
the galaxy type and star-formation history. iBendo et al.l 
12008) find that the 8/160 /im ratio in 15 nearby face-on 
spirals suggests that the PAHs are associated with the 
diffuse cold dust that pro duces most of the 160 /im emis- 
sion. On the other hand, iHaas et all ()2OO20 find that the 
8 ^m feature, across a range of galaxy type and current 
star formation rate, is associated with peaks of 850 /xm 
surface brightness, which originate in molecular regions. 
The distribution of PAHs in a galaxy is one parameter 
that will help determine what fraction of the PAH lumi- 
nosity aries from the reprocessing of UV photons from 
young, hot stars versus the general galactic distribution 
of B stars. This distinction is crucial in usi ng PAH emis- 
sion a s a tracer of current star-formation (iPeeters et al.l 
12004 iCalzetti et"ani2007D . To further understand the 
distribution of PAHs in the SMC, we explore the corre- 
lation of PAH fraction with 160 emission and the dust 
mass in Figure [HI In this Figure, we see that the PAH 
fraction is correlated with dust surface density (Md) and 
160 /im emission, but only weakly correlated with H I 
column (note that the H I column is shown in a linear 
scale while M^ and 160 /im emission are shown on a 
logarithmic scale). 
The correlation of the PAH fraction with dust surface 



density but not H I reffects the fact that PAHs are not 
uniformly distributed in the SMC. Regions with PAH 
fraction greater than 1% in general have M^ > 10^ 
Mq kpc~^. However, regions with dust surface densities 
above this level also tend to contain molecular gas, so 
the dust surface density and H I column no longer track 
each other because of the presence of H2 (jLerov et al.l 
l20f)l . For this reason, we see at best a weak correlation 
of gpAH with neutral hydrogen column, but a stronger as- 
sociation with CO emission. Bolatto et al 2010 (in prep) 
have used the MIPS observations of the SMC from S^MC 
and the SAGE-SMC survey to map the distribution of 
molecular gas as inferred from regions with "excess" dust 
surface density relative to the co lumn of neu t ral hy dro- 
gen, us ing the same techniques as lLerov et al.l (|2007[ ) and 
iLerov et al. (2009). In Figure [TSl we show their map of 
the H2 column density overlayed with a contour at (7pah= 
1% and a contour of 3cr CO emission. Although we must 
use caution in comparing the detailed distribution of CO 
to gpAH since the CO is at ~ 4 times lower resolution, 
it is interesting to note that there are regions with high 
molecular gas columns without CO and with low PAH 
fractions, for instance in the northern region of the SW 
Bar. This may indicate that the radiation field in these 
regions is affecting both CO and PAHs. 
We have so far shown that PAH fraction is high in re- 



12 



Sandstrom et al. 





Fic;. 4.- - Histograms illustrating the variation of the radiation field paramters for the photofit models (in black) and photospectrofit 
models (in gray). The mean radition field (U) is shown in the lower right panel. Despite the redistribution of the radiation field, the mean 
field is essentially unchanged in the two best fit models. 




Fig. 5. — A comparison of fpoR for the photofit and photospectrofit models. The dashed line on the histograms illustrates the mean of 
the fpDR values. On average fpoR increases by 0.074 when the spectroscopic information is included in the fit. 
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Fig. 6. — Comparison of some representative SEDs and photofit models. The positions of these SEDs are marked with square symbols 
on Figure [7] The gray line appears in each panel for comparison and shows one of the highest gpAH values from the model fits from the 
SW Bar which is marked on Figure [7] with a black square. The location of the SEDs are listed in the top left of the plot. The measured 
photometric points are shown with error bars and the synthetic photometry for the best fit model is shown with a filled circle. These panels 
illustrate the range of qpAH values we see in the SMC. 




Fig. 7. — This map shows the gpAH values from our fit to the photometry at each pixel in the mapped region (40" resolution). The outer 
white boundary shows the overlappi ng coverag e of the MIPS and IRAC mosa ics. The white asterisks show the locations of the stars with 
measured UV extinction curves (Go rdon et al.l[20 031 discussed in Section 14.51 The white squares show the locations of the SEDs plotted 
in in black in Figure [6] The black square in the SW Bar shows the location of the SED plotted in gray in each panel of Figure [6] 



gions of active star-formation, associated with the pres- 
ence of CO and molecular gas. PAHs can also be de- 
stroyed in regions of active star-formation by the intense 
UV fields produced by massive stars o r in the H 11 region s 
themselves by chemistry with 11+ (jGiard et al.l 1199^ . 
Figure [m shows the MCELS map of Ha in the SMC 
overlayed with the 1% contour of gpAH- From this com- 
parison it is clear that the regions of high PAH fraction 
are typically on the outskirts of H 11 regions (i.e. the high 
9PAH regions and the H 11 regions are not co-spatial). In 
particular, the region around N 66, the largest H 11 region 
in the SMC, has a very low PAH fraction. We will discuss 
the influence of H 11 regions and massive star-formation 
on the PAH fraction further in Section [531 

4.4. The PAH Fraction in SMC Bl # 1 

The molecular cloud SMC Bl 1 was the first location 
in the SMC where the emission from PAHs was identified 
(|Reach et al.|[200l . This region has been studied exten- 
siv ely and the PAH em ission spectrum has been modeled 
bv lLi fc Draind H2002f ) who found that the PAH fraction 
in SMC Bl (gpAH~ 1.6%) was 8 times higher than the 
average frac ti on in the Bar (qpah"^ . 2%). In addition, 
IReach et all (|2000D and iLi fc Draint^ (|2002D noted un- 
usual band ratios of the 6.2, 7.7, and 11.3 /zm features. 
In Figure [T7] we show the best fit models for the spectrum 
and photometry of SMC Bl. We fin d (/pah-- 1.2 d = 0.1% 
in relatively good agreement with iLi fc Draind (|2002j ) 



considering the differences in resolution between our re- 
spective datasets. We also reproduce the distinctive band 
ratios (11.3 and 6.2 features high compared to the 7.7 fea- 
ture) seen with ISO. SMC Bl does not have a uniquely 
large PAH fraction compared to other locations in the 
SW Bar. 

4.5. PAH Fraction and the 2175 A Bump 

On Figures [7] and [TH] we show the locations of the five 
stars with measurements of their UV extinction curves 
with asterisks. Of these stars, only one shows a 2175 A 
bump in its extinction curve: AzV 456, which unfortu- 
nately lies just outside the boundaries of the map. The 
lack of a bump in the remaining stars has been inter- 
preted as evidence for a low PAH fraction along those 
lines-of-sight, assuming that PAHs are the carrier of the 
bump (Li &: Drainc 2002). To test this assertion, we list 
in Table [6] the measured values of E(B— V) for the five 
stars with extinction curves from [Gordon et all (pOOl 
along with the E(B— V) and gpAH we calculate from 
our photofit model results at those positions. For the 
[Drainc & Li (2007) models, the dust mass surface den- 
sity yio and E(B— V) are related by a constant value of 
2.16 X 10-6 mag (M© kpc'^)"!. 

The comparison of the E(B— V) values shows that the 
stars are indeed behind the majority of the dust along 
those lines of sight, and that the gpAH values are slightly 
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Fig. 8. — A map of 8 /im emission overlayed with the 1% contour of gpAH- There are a number of regions that are very bright at 8 //m 
that do not have high PAH fractions, particularly N 66 and the northern part of the SW Bar. 




qpAH (%) 

Fig. 9. — This figure shows a two-dimensional histogram of the 8/24 ratio vs gpAH overlayed with the binned average of the 8/24 ratio 
and error bars representing the standard deviation of the scatter at ea<;h gpAH bin. The color scale shows the number of points in each bin. 



16 



Sandstrom et al. 



Logio(Ub^,) 



0.0 0.5 1.0 1.5 0.0 0.1 0.2 0.3 




Fig. 10. — Maps of the average radiation field (U) and PDR fraction {/pdr) from the photofit models. We have overlayed one represen- 
tative contour of the MCELS Ha image to illustrate the locations of the brightest H II regions in the SMC. The correspondence between 
fpDR and the location of H II regions is very good, as expected. 




LogioCMo) (Mo/pc^) Leroy et al. (2007) 

Fig. 11. — A tw o-dimensional histogram showing the comparison of the dust mass surface density M-d from our study with that of 
ILerov et al.l 1120071). The gr ay-scale, with a linear stretch, shows the density of points in the plot. Our Mjj is approximately 70% of that 
found bv ILerov etaLl 120071 ). using the same MIPS data but different methodolgy. We note that the scatter at low surface dens ities likely 
relates to the fact that we performed our analysis at 40" resolution but convolved to 2.6' for comparison with ILerov et ahl 112007 ) and we 
would expect higher signal-to-noise if the analysis had been performed in the opposite order. 
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Fig. 12. — In this figure we show the 5pah map from Figure [7] with contours of 3(t CO (J=l— 0) emission from the NANTEN survey 
overlayed. The thin black line represents the coverage of the NANTEN survey. The NANTEN observations have a resolution of 2.6'. 



above the SMC average of 0.6% (see Section l4?2|) . A more 
detailed analysis of these lines-of-sight will be presented 
in a follow-up paper with targeted IRS spectroscopy to 
study the PAH emission in these regions. Since we do 
find that regions of high <7pah tend to be associated with 
molecular gas, it may be the case that assuming the dust 
and PAHs are uniformly mixed along each line of sight 
does not hold. In that case, the comparison of E(B— V) 
values may not be a good indicator of whether these stars 
should show the 2175 A bump in their extinction curves. 
In addition, our angular resolution is not high enough 
in these maps to directly observe the structure of the 
dust emission in the vicinity of these stars, so we can 
not definitively test the PAII-2175 bump connection. We 
note that one of the stars lies near N 66 in a region with 
very low ambient PAH fraction. Some of the other stars 
may fall in voids in the PAH distribution, but higher 
angular resolution is necessary to understand the line of 
sight towards these stars. 

5. WHAT GOVERNS THE PAH FRACTION IN THE 

SMC? 

iDraine et al.l (|2007[ ) determined the PAH fraction in 



the SINGs galaxy sample using identical techniques to 
what we have done here. They observed a e f fect v ery 
similar to what was seen bv lEngelbracht et al.l (|2005[ ) in 
that at a metallicity of 12+ log(0/H) ~ 8.0 there is an 
abrupt change in the typical PAH fraction or 8/24 
ratio. The SMC lies right at this transition metallicity, 
so we hope to gain some insight into the processes at 
work by studying its PAH fraction. There is still a great 
deal of uncertainty as to how PAHs form and how they 
are destroyed, let alone how those processes balance in 
the ISM. In the following sections we address some of 
the proposed aspects of the PAH life-cycle and elucidate 
what we can learn about them from the SMC. 

5.1. Formation by Carbon- Rich Evolved Stars 

The formation of PAHs in the atmospheres of evolved 
stars is a well established hypothesis for the source of 
interstellar PAHs. PAH emission bands have been ob- 
served in the spectrum of the carbon-rich post-AGB stars 
where the radiation field increases in hardness and in- 
tensit y, more effectively exciting the PAH emission fea- 
tures (|Justtanont et al.l I1996D . Despite these observa- 
tions of PAH formation in carbon-rich stars, it remains 
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Fig. 13. — Histogram of gpAH in lines-of-sight with and without detected CO emission from the NANTEN survey. 




Fig. 14. — Two-dimensional histograms of the dust mass surface density, 160 /im surface brightness, neutral hydrogen column and CO 
integrated intensity as a function of gpAH overlayed with the binned average. The error bars show the standard deviation of the scatter in 
each bin of gpAH- The gray scale represents the logarithm of the number of points in each bin. 
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Fig. 15. — A map of molecular gas column density inferred from excess dust emission at 160 fira from Bolatto et al 2009 (in prep) 
overlayed with the 3cr contour of CO emission in green and the 1% contour of (/pah in red. The coverage of the CO survey is shown as a 
thin gray line. 



TABLE 6 

Comparison with Extinction Curve Measurements 



Star 


R.A. 


Dec. 


E(B-V) (Gordon et al. 2003) 


E(B-V) (This Study) 


9PAH 




(J2000) 


(J2000) 


(mag) 


(mag) 


(%) 


AzV 18 


0'^47™12= 


-73°06'33" 


0.167±0.013 


0.25±0.07 


0.4± 0.1 


AzV 23 


0'^47™39= 


-73°22'53" 


0.182±0.006 


0.30±0.10 


0.8± 0.1 


AzV 214 


0h58'"55= 


-72°13'17" 


0.147±0.012 


0.34±0.13 


0.7± 0.2 


AzV 398 


1^06^10= 


-71°56'01" 


0.218±0.024 


0.33±0.09 


0.8± 0.1 


AzV 456 


1^10™ 56= 


-72°42'56" 


0.263±0.016 







Note. — 



The E(B— V) calculations arc described in Scction |4.5l 
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Fig. 16. — MCELS Ha image with the 1% contour of gpAH overlayed. The regions of high PAH fraction are typically on the outskirts of 
H H regions. 



to be shown that they inject PAHs into the ISM at the 
level needed to explain the abundance observed. This, 
of cou rse, is a gen eral problem in the "stardust" sce- 
nario (Iprainell2009f). Recentl y, in the Large Magellanic 
Cloud. TMatsuura et al.l ()2009f ) have performed a detailed 
accounting of the dust enrichment of the ISM by AGB 
stars and find a significant deficit compared to the ob- 
served ISM dust mass. 

Assuming that th e " stardu st" racture is correct, 

iGalliano et aH (|2008[ ) and lDwekl (|1998[ ) hypothesize that 
the low fraction of PAHs in low metallicity galaxies re- 
flects the delay (~ 500 Myr) in the production of carbon 
dust from AGB stars relative to silicate dust from core- 
collapse supernova. For the SMC in particular, this pic- 
ture has a number of issues. Most importantly, there is 
evidence that the SMC formed a large fraction of its stars 
more than 8 Gyr ago followed by a subsequent burst of 
star formation 3 Gyr ago (Harris & Zaritsky 2004). This 
long time scale makes the delayed PAH injection into the 
ISM by AGB stars an unlikely explanation fo r the current 
observ ed PAH deficiency. Recent work by iSloan et"all 
(j2009( ). for example, argues that carbon- rich AGB stars 
in low metallicity galaxies can start contributing dust to 



the ISM in - 300 Myr. 

To evaluate the relationship between the current distri- 
bution of PAHs and the input of PAHs from AGB stars, 
we show in Figure[T9]the PAH fracion overlayed with con- 
tours of the density of carbon AGB stars. These contours 
are created using the 2MASS 6X point source catalog 
towards the SMC, selecting carbon st ars using the tech- 
nique described in iCioni et alJ (|2006l ). The distribution 
of carbon stars follows the observed "sphe roidal" popula- 
tion o f older sta rs in the SMC very well (jZaritskv et alJ 
l2000t ICioni et "a l. 2000). The (?pah map, however, has 
no clear relationship to the distribution of carbon stars. 
This is perhaps to be expected since the distribution of 
ISM mass does not follow the stellar component either. 
However, we note that a different conclusion regarding 
the PAH fracion compared t o the distribut i on of AGB 
stars was recently reached by 'Parad is et al.l ()2009D , who 
modeled dust SEDs acro ss the Large Magel lanic Cloud 
using the SED models of lDesert etaJ] (|1990[ ). They find 
evidence that the fracion of PAHs in the LMC is highest 
in the region of the stellar bar, where the concentration 
of AGB stars peaks. This could be an accidental coinci- 
dence, or it may refiect methodological differences or it 
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Fig. 17. — Best fit models to the photometry and spectroscopy in SMC Bl # 1. The gray points show the IRS spectrum of SMC Bl, 
while the diamonds show the photometry from IRAC and MIPS. 



could indicate a different dominant mechanism of PAH 
formation in the SMC and LMC or more efficient de- 
struction of PAHs in the SMC. 

Although the distribution of PAHs may not resemble 
that of AGB stars, we would expect the PAHs injected 
by those stars to be present in the diffuse ISM and not 
preferentially in molecular clouds. As we have shown, 
however, the diffuse ISM PAH content is very low in 
the SMC. As such, in order to reconcile the pathway 
of PAHs from AGB stars to the diffuse ISM to molec- 
ular clouds, we must hypothesize a recent event, occur- 
ing after the condensation of the current generation of 
molecular clouds, which essentially cleared the diffuse 
ISM of AGB-produced PAHs while leaving the shielded 
regions of molecular gas with high PAH fractions. Ob- 
servations of giant molecular clouds (GMCs) in the Mag- 
ellanic Clouds and other ne arby galaxies sugg est that the 
CMC hfetime is ~ 25 Myr (|Fukui et al.lll999HBlitz et al.l 
I2007D . so the event in question would have had to occur 
within the last 25 Myr or so. In the absence of such an 
event (which we will investigate further in a Section [531) , 
our observation of low PAH fraction in the diffuse ISM, 
higher PAH fraction in molecular clouds, and the lack of 
relation between the PAH fraction and the distribution 
of AGB stars is strong evidence against AGB stars being 
the dominant force behind the fraction of PAHs in the 
SMC. 

5.2. Formation and Destruction of PAHs in Shocks 
and Turbulence 

The shocks created by supernova explosions have a dra- 
matic effect on the content and size distribution of dust 
in the interstellar medium. Upon encountering a shock 
wave, grains can be shattered via collisions or sputtered 
by hot gas behind the shock or by motion of the grain 
through the post-shock medium. Because of their small 
mass-to-area ratio, PAHs are well coupled to the gas and 
do not acquire large relative velocities after the passage 
of a shock, so they will primarily be sputtered only in hot 



post-shock g as behind fast (v > 200 km s~^) shocks. Cal- 
culations by iJones et"aLl ()1996[ ) suggest that grain shat- 
tering could in fact be a net source of PAH material for 
shocks between 50 and 200 km s~^, converting ~ 10% 
of the initial grain mass into small PAH sized fragments. 
Thus, it is not immediately obvious what the net effect 
of interstellar shocks on the fraction of PAHs will be. 

Some studies attribute the low PAH fraction in low 
metallicity galaxie s to efficient dest r uction of PAHs by 
supernova shocks. lO'Halloran et all (|2006f) found a cor- 
relation between decreasing PAH emission and increased 
supernova activity as traced by the ratio of mid-IR lines 
of [Fc II] and [Ne II] . Two issues with this interpretation 
are that the mid-IR lines are tracing current supernova 
activity, which only affects the PAH fraction in the im- 
mediate vicinity of those remnants, and an increased su- 
pernova rate will have recently been related to a higher 
UV field produced by the massive stars, so it is difficult 
to disentangle the effects of shocks versus intense UV 
fields. 

The distinctive nature of the SMC extinction curve 
might point toward the influence of supernova shocks on 
the dust grain size distribution. Four of the five extinc- 
tion curves show similar characteristics: a lack of the 
2175 A bump (the carrier of which is most likely PAHs 
iLi fc Drain jboOlf ) and a steeper far-UV rise indicating 
more small dust grains relative to the Milky Way extinc- 
tion curve (jCordon et all 120031: iCartledge eFall I2005D . 
Magellanic-type extinction has been seen to arise in 
the Milky Way as well along the sightline towards HD 
204827, which may be embedded in dus t associated with 
a supernova shock (jValencic et al.|[2003[ ). However, there 
are many other viable interpretations to explain the pro- 
portion of small dust grains, including inhibition of grain 
growth in dense clouds. 

In Figure [20] we show a map of the H I v elocit y 
dispersion in the SMC from iStanimirovic et al.l ([1999). 
The velocity dispersion here mostly traces the regions 
where there are more than one velocity component along 
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Fig. 18. — The dust mass surface density and E(B— V) values derived from our photofit models at each pixel overlayed with contours of 
neutral hydrogen column density from Stanimirovic ct al. (199$|) at 6, 10, 14 and 18 X lO^'^ cm~^. The locations of the five stars in the 
SMC with extinction curves from lGoTdon et al.i t2003i) are shown with asterisks. Table [6l shows a comparison of the E(B— V) measured for 
those stars with the total line-of-sight E(B— V) we calculate from the photofit model results. 



the li ne of sight, particularly hiKhlighting the regions 
where [stanimirovic et al.l ()1999f ) find evidence for super- 
giant shells in the H I distribution. We show the ap- 
proximate locations of two of their shells that overlap 
our map. If supernovae are the source of these shells, 
iStanimirovic et al.l ()1999[ ) find that - 1000 supernovae 
are required to account for the kinetic energy and the 
shells have dynamical ages of ^ 20 Myr. The supergiant 
shells provide indirect evidence for the effects of super- 
novae on the ISM. On Figure [50] we also mark with green 
crosses the locations of young 1000 — 10000 yr) su- 
pernova remnants identified in the Australia Telescope 
Compact Array survey of the SMC (Payne_et al. 2004 ). 

There is no clear trend relating (/pah to the boundaries 
of the supergiant shells, although this may be an effect 
of depth along the the line of sight. The middle region of 
the Bar and Wing, which is essentially devoid of PAHs 
is covered by one of the shells, but the SW Bar, which 
hosts the largest concentration of PAHs in the SMC is 
covered as well. Although there is an anticorrelation be- 
tween the young SNRs and large PAH fraction, the dis- 



tribution of remnants closely follows that of the massive 
star- forming regions (compare Figures 1161 and [20 ]) . For 
this reason, it is difficult to draw strong conclusions as to 
whether the supernova shocks or the UV fields and H II 
regions created by their progenitor stars were responsible 
for destroying PAHs in these areas. 

Similar to shocks, regions of high turbulent velocity 
may alter the size distribution of du st grains via shatter- 
ing. IMiville-Deschenes et al.l (|2002[ ) studied a region of 
high latitude cirrus and argued that variations in small 
dust grain and PAH fractions coul d be related to the tur- 
bulent velocity field in the region. iBurkhart et"all pOlol ) 
have studied turbulenc e in the ISM of the SMC using 
the H I observations of IStanimirovic et al.l ()1999l ). They 
present a map showing an estimate of the sonic Mach 
number of turbulence in the SMC based on higher order 
moments of H I column density and the results of nu- 
merical simulations. This map is quite distinct from the 
velocity dispersion map shown in Figure [20] which mostly 
highlights the presence of bulk velocity m otions along the 
line of sight. In Figure I2l1 we show the IBurkhart et aD 
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Fig. 19. — Map of gpAH ovcrlaycd with contours of carbon AGB star density determined from the 2MASS 6X point source catalog using 
the selection criteria of Cioni ct al. (2006). The contours are labeled with the density of carbon stars per square degree. There is no obvious 
correspondence between the distribution of carbon stars and the fraction of PAHs. 



(|201Clf) Mach number map overlayed with the contours 
of (?PAH- We find that there is no strong correlation be- 
tween tubulent Mach number and (Jpah, although in gen- 
eral the regions of high (^pah are near depressions in the 
Mach number distribution, opposite what we would ex- 
pect if grain shattering in turbulence was a major source 
of PAHs. However, this simple comparison does not pro- 
vide information about the amount of dust affected by 
different levels of turbulence (e.g. some of the low Mach 
number regions may represent a very small fraction of 
the total ISM). To account for this fact, we also calculate 
the dust mass surface density weighted Mach number in 
regions with g'pAH> 1% and < 1% to be 0.8 ± 0.3 and 
1.0 ± 0.4, respectively. Although there is no clear asso- 
ciation between high sonic Mach number and gpAH from 
this comparison, it is worth noting that if PAHs are pri- 
marily associated with molecular gas, the Mach number 
derived from H 1 observations may have little relevance 
to the creation or destruction of PAHs in the SMC. 

In general, we do not see clear cut evidence that turbu- 
lence or shocks are the major drivers of the PAH fraction 
in the SMC. Young SNRs are found around H 11 regions. 



making it difficult to separate the effects of radiation field 
and shocks in the destruction of PAHs in those regions. 
The supergiant H 1 shells observed in the SMC may be a 
tracer of where supernova shocks have seriously affected 
the ISM, but we see peaks of gpAH within their bound- 
aries barring a line-of-sight depth effect. 

5.3. Destruction of PAHs by UV Fields 

In low metallicity galaxies, the dust to gas ratio is de- 
creased (Lisenfeld & Ferrara 1998) and the effects of the 
UV field from regions of massive star formation can be 
spread over a much larger area. In addition, the de- 
creased metallicity may lead to harder UV fields because 
of the lower line blanketing in stellar atmospheres. These 
changes can be t raced by the ra tios of mid-lR emission 
lines and work bv lMaddenI (|2000D and others have shown 
that ra diation fields a r e hard er in low metallicity environ- 
ments. IGordon et al.l (j2008[ ) studied star-forming regions 
in M 101, over a range of excitation conditions and metal- 
licities. They found that the strength of PAH emission 
decreases with increasing ionization parameter, suggest- 
ing that processing by the UV field from these regions 
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Fig. 20. — SMC H I velocity dispersion from'Stanimirov ic et al.] I I1999I) . The approximate locations of the supergiant H I shells identified 
by Stanirniroyic et al. ( 1999) are shown with a dashed white line. The velocity dispersion mainly traces bulk motions of gas along the line 
of sight, particularly highlighting the regions of these shells. We show the locations of young supernova remnants identiiied in the ATCA 
survey of the SMC with black crosses. 



was the major force behind the changing PAH fraction. 

The exact mechanism of PAH destruction by UV fields 
is not entirely clear. Small PAHs can be destroyed 
by the ejection of an acetylene gro up upon absorbing 
a UV photon (jAllain et al.l Il996al) , but PAHs larger 
than ^ 50 carbon atoms are relatively stable to these 
effects under a range of UV field conditions. If the 
PAHs are partially dehydrogentated or hi ghly ionized, 
they a re more susceptible to destruction (jAllain et alj 
ll996bD . so one possibility is that PAHs in low metallicity 
galaxies tend to be more highly ionized, dehydrogenated 
or smaller than their counterparts in higher metallicity 
galaxies. However, there is little evidenc e from the spec- 
tra of low ineta llicity galaxies (Engclbra cht et all 120081 : 
ISmith et a"l][2007i ) that PAHs are different at low metal- 
licity. 

There is also a distinction to be made between the 
destruction of PAHs by intense radiation fields in the 
immediate vicinity of H H regions and a global decrease 
of the PAH fraction in the galaxy. We see evidence for 
PAH destruction near H II regions in our map (see Fig- 



ure [16]). But can UV fields from massive star forming 
regions be responsible for the low gpAH over the entire 
galaxy? There is evidence that star-forming regions in 
irregular galaxies can "leak" a large fraction of their ion- 
izing photons. Observations of emission line ratios from 
the diffuse ionized gas (DIG) in star-forming dwarfs in- 
dicate that dilution of radiation from a central massive 
star-forming region which leaks a significant fra ction of 
its io nizing photons is a likely source for the DIG (jMartinl 
[1993). In the SMC, models of N 66 suggest that ~ 45% 
of the ionizing photons esc ape the H II region and go on 
to ionize the diffuse ISM (jRelaiio et all 120021 ). Thus, it 
is at least plausible that UV fields may be an important 
driver of PAH destruction over the entire galaxy. 

To evaluate the possibility that the low gpAH in the 
SMC is due to destruction by UV f ields, we examine 
the re solved star- formation history of [Harris fc Zaritskvl 
(|2004f ) to search for recent star-forming events that could 
have affected the PAH fraction in the diffuse ISM through 
their UV fields. In Figure [22] we sh ow five panels illus- 
trating the [Harris fc Zaritskvl ([20041 ) results regridded to 
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Fig. 21. — Map of estimated turbulent Mach number from Burkhart et al (2009) in prep, overlayed with the 1% contour of (/pah- The 
Mach number map has an angular resolution of 30'. 



match our map of qpah- These panels show the total 
star formation rate for bins older than 1 Gyr, from 400 
Myr to 1 Gyr ago, from 100 to 400 Myr ago, from 25 to 
100 Myr ago and younger than 25 Myr. 

In general, panels which show the star-formation rate 
in bins older than 100 Myr tend to follow the spheroidal 
distribution of older stars, as traced by the carbon stars 
shown in Figure \W[ The 25 — 100 Myr panel shows that 
star formation occurred along the Bar, coinciding with 
the region in our gpAH map that is devoid of PAHs. The 
most recent star formation, which would overlap the time 
when the current generation of molecular clouds con- 
densed, seems to be mainly associated with the outer- 
most edges of the Bar, near where we see PAHs. 

PAH destruction as a side effect of massive star for- 
mation is difficult to disentangle since the UV fields and 
subsequent supernovae only affect the dust over a short 
period of time surrounding the star formation event. Be- 
tween 25 and 100 Myr ago in the SMC, the star-formation 
along the Bar, which most likely relates to the supergiant 
bubble seen in the H I overlapping that location, may 
have cleared the region of PAHs by some combination 
of UV fields and shocks. However, most of this activity 



likely occured before the condensation of the molecular 
cloud s (~ 25 Myr > ago: iFukui et aLl ll999l: IBlitz et al.l 
I2007t IKawamura et al.ll2009l )" that currently have a high 
PAH fraction, so there is still difficulty in reconciling the 
high gpAH in molecular clouds and the low (Jpah in the 
diffuse ISM by destruction alone. 

5.4. The Formation of PAHs in Molecular Clouds 

One of our primary observations in the SMC is that 
high PAH fractions occur along lines of sight through 
molecular gas. We argue that such a situation could 
arise in two scenarios: 1) AGB stars enrich the diffuse 
ISM with PAHs, part of which is then incorporated into 
molecular clouds. A subsequent event (e.g. SN) clears 
the PAHs in the diffuse ISM. Or, 2) PAH formation oc- 
curs in the molecular cl ouds. These t wo sc enarios are 
not mutually exclusive. 'Paradis et alj ()2009( ) found en- 
hanced PAH fraction in the LMC both in the stellar bar, 
which hosts the highest concentration of AGB stars, and 
in molecular clouds. It is possible that in the diffuse ISM 
of the SMC, AGB produced PAHs are rapidly destroyed 
and all we observe are PAHs that formed in molecular 
clouds. 





Fig. 22. — Resolved star-formation history of the SMC from lHarris &: Zaritskvl 120041) . The time-ranges are hsted in the upper left corners 
of each image. 
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iGreenberg eFal] (|2000l ) propose a scenario by which 
PAHs could be formed in dense clouds. They argue that 
a layer of ices and organic material forms on grains in 
dense clouds and the photoprocessing which occurs in the 
transition from the dense cloud back to the diffuse ISM 
forms PAHs, a phenomenon they explore through labo- 
ratory experimentation. There is also observational evi- 
dence that hints at some of the p r ocessi ng P AHs undergo 
in den se clouds. iRapacioli et al.l ()2005[) and iBerne et al.l 
(|2007| ) show variations in the PAH spectrum in photo- 
dissociation regions consistent with PAHs being in clus- 
ters or embedded in a carbonaceous matrix in the sur- 
rounding molecular gas and emerging as free- flying PAHs 
closer to the exciting star. 

While it would not be out of question that a burst of 
star-formation ^ 25 Myr ago could have a dramatic effect 
on the global diffuse ISM of the SMC, a more likely sce- 
nario for the concentration of PAHs in molecular clouds 
is that they form there. Assuming that is the case there 
are a number of interesting observations that we can 
make about the formation of PAHs in molecular clouds 
in the SMC. First, the correspondence between the CO 
and gpAH is very good, better even then the correspon- 
dence between the PAHs and the Bolatto et al 2009 (in 
prep) map of molecular hydrogen inferred from excess 
dust emission. This suggests that the conditions for the 
formation of CO and the formation of PAHs are simi- 
lar. It may be that the UV field in the outskirts of the 
clouds, where H2 can self-shield but CO is dissociated, 
can disrupt the growth of mantles on the grains, prevent- 
ing a process that creates PAHs. Alternatively, because 
the growth of mantles requires relatively high densities, 
the presence of CO may allow the gas to cool and be- 
come denser, increasing the rate of gas-grain reactions. 
This scenario could be tested by comparing the distribu- 
tion of CO and PAHs in a sample of irregular galaxies. 
iLerov et all (|2b09^ found that CO in the N 83 region 
of the SMC is found where Ay > 2 through the cloud 
(Ay > 1 at cloud center). If CO survives and PAHs form 
at a similar extinction, the decrease in the PAH fraction 
and the deficit of CO emission at low metallicity may be 
related to the scarcity of regions where this condition is 
met. 

Second, we observe that the typical gpAH in SMC 
molecular clouds is ~ 1 — 2%. This is lower than the 
average Milky Way value (which is for the diffuse ISM) 
and may represent some limit on the efficiency of creating 
PAHs in molecular gas in the SMC. In the low metallic- 
ity galaxies studied bv iDraine et al.l (|2007[ ). those with 
metallicities below 12-1- log(0/H) ~ 8 had a median PAH 
qpAn of 1% with large scatter. This scatter could be due 
to the filling factor of molecular gas. If PAHs are forming 
primarily in molecular clouds, the differences observed in 
the PAH fraction in low versus high metallicity galaxies 
may be related to the ability of PAHs to survive and ac- 
cumulate in the diffuse ISM and/or to the efficiency of 
PAH formation in molecular gas, which can be decreased 
by the lower carbon abundance and the more pervasive 
UV fields in low metallicity environments. 

Finally, in the scenario where PAHs are destroyed in 
the diffuse ISM and formed in molecular clouds and 
AGB stars, the abrupt transition at a metallicity of 12-|- 
log(0/H) ~ 8 would represent mostly a change in the ef- 
ficiency of PAH destruction in the diffuse ISM. If PAHs 



are destroyed efficiently in the diffuse ISM of these sys- 
tems, all that is left are the PAHs created in molecu- 
lar clouds. A resolved study of the PAH abundance in 
the SINGs galaxies using higher resolution far-IR and 
millimeter-wave CO observations from Herschel may be 
able to show whether the PAH fraction in the diffuse ISM 
versus molecular clouds is changing across the transition 
metallicity. 

6. SUMMARY & CONCLUSIONS 

We present results of fitting the IDraine fcTil ()2007l ) 
dust models to SEDs and spectra obtained from S^MC 
and S'^MC. Our major results are as follows: 

1. Comparisons of best-fit dust models in regions with 
overlapping photometry and spectroscopy demon- 
strate that the IRAC and MIPS SED estimator for 
9PAH does not appear to be appreciably biased de- 
spite the absence of information between 8 — 24 and 
24 — 70 /.tm. When 5 — 38 jum IRS spectroscopy is 
added to the fitting constraints, the resulting "pho- 
tospectrofit" models yield (/pah values that are only 
slightly larger than the "photofit" model estimates. 
The photospectrofit models have more dust at in- 
termediate (T w 60K) temperatures, and less dust 
with T w 200 K. The estimate for the average radi- 
ation intensity scale factor U and the dust surface 
density M^ are nearly unaffected. 

2. The PAH fraction in the SMC is low and variable. 
As a fraction of the total dust mass, the highest 
PAH fractions we observe are about half of the 
Milky Way value, but most of the galaxy has PAH 
fractions at the lower limit of our models {qpAH— 
0.4%), an order of magnitude lower than the Milky 
Way value. The average (qpAn) ^ 0.6% in the 
SMC. This is consistent with the earli er estimate of 
a ver y low PAH fraction in the SMC ()Li fc Draind 
120021 ) based on IRAS and DIRBE photometry. 

3. The 8 /im emission alone does not trace the (Jpah 
weh in the SMC. We find that the 8/24 ratio is cor- 
related with (jpAH and agrees with what one would 
predict from the metallicity of the SMC an d the ob- 
served trends in lEngelbracht et"all (|2005 h. but for 
a given gpAH there is a wide range of 8/24 ratios, 
depending on the intensity of the local starlight 
heating the grains. 

4. The metallicity of the SMC places it at the transi- 
tion where lower metallicity galaxies show a de- 
ficiency of PAHs and higher metallicity galaxies 
show approximately Milky Way level PAH frac- 
tions. If the SMC is typical of galaxies in this re- 
gion, the transition seems to represent a decrease in 
the filling factor of PAH-rich regions rather than a 
uniform decrease in the PAH fraction throughout 
the galaxy, although even the highest PAH frac- 
tions are still well below MW levels. 

5. Lines-of-sight with molecular gas have an average 
qpAH of ^ 1%, while lines of sight through the dif- 
fuse ISM have qpah at least a factor of two less. 
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6. We evaluate the various proposed drivers of the 
deficit of PAHs at low metallicity. The distribu- 
tion of PAHs in the SMC does not follow the carbon 
AGB star distribution, the regions of high turbu- 
lent Mach number, or the location of supergiant 
shells in the ISM. 

7. The low PAH fraction in the diffuse ISM versus the 
high PAH fraction in molecular regions leads us to 
propose that PAHs may be forming in molecular 
clouds and/or a recent event (perhaps related to 
star formation events in the last ~ 25 Myr) de- 
stroyed a large fraction of the PAHs in the diffuse 
ISM. 

8. We surmise that the global PAH fraction at low 
metallicities is a reflection of the amount of gas 
in these systems that is found at high extinction 
(Ay > 1). 
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APPENDIX 
Cross- Calibration of MIPS, IRAC and IRS 

Because we combine photometry and spectroscopy in our model fitting, we investigated whether there were systematic 
differences between the calibrations of the IRS and MIPS/IRAC instruments. The SSC finds the cross-calibration agrees 
to within 10% for the three instruments. However, most of these tests have been carried out on stellar sources, whereas 
our observations are primarily concerned with extended emission. We have used the extended source corrections for 
IRAC determined bv iReach et al.l ()2005| ) and the slit loss correction function for IRS distributed with Cubism. 

We have constructed 8 and 24 jjwi maps from t he spectral cubes using the most recent versions of the spectral 
response curves available from the SSC (for IRAC; iHora et all 120081) . We then compare the IRS 8 and 24 /xm maps 
to the same regions in the IRAC and MIPS mosaic after smoothing to eliminate any small PSF differences between 
the two instruments. Calibration differences between the instruments present themselves as a non-unity slope in the 
comparison of IRS 24 /im to MIPS 24 /im and IRS 8 /im to IRAC 8 /xm. For the 8 /im comparison, we find that the 
IRS photometry is higher than the IRAC photometry by factor of 1.1, on average, although the slope in individual 
regions varies between 1.04 to 1.14. At 24 /im the IRS photometry is lower than the MIPS photometry by a factor of 
0.94 on average, varying between 0.91 and 0.96 region to region. Individually, these differences are within the quoted 
cross-calibration limits stated by the SSC, however without correction they introduce scatter in the comparison of 
spectroscopically and photometrically determined (/pah- 

In addition to the slope, there is an additive offset between the IRS and MIPS/IRAC photometry that is too large 
to be an issue with the MIPS/IRAC mosaic foreground subtraction described previously (additive offsets are ~ 0.1 — 1 
MJy sr-i at 24 //m and ~ 0.01 - 0.1 MJy sr'^ at 8 urn). The source of this additive offset is not definitively known, 
but we consider it very likely that it results from slight under or over subtraction of the zodiacal light foreground as 
determined by SPOT for the IRS cubes. The uncertainties of the DIRBE model are on the order of a few percent at 
the DIRBE wavelengths, and particularly at 24 /im where the zodiacal foreground is very bright (on the order of tens 
of MJy sr~^) and a few percent error easily translates into offsets on the order of a few tenths of a MJy sr~^. In order 
to eliminate these offsets, we must make some assumption about the spectral shape of the correction. To do this we 
fix the "off" position zodiacal spectrum and determine the percentage by which the "map" position spectrum must be 
off to produce the offsets we observe. We do this for each AOR in the dataset and then add in the resulting correction 
spectrum. The correction is less than 5% in all cases. 

The multiplicative differences in the calibration between IRS and MIPS most likely result from the extended source 
calibration in IRS, since MIPS does not require such a correction. For the IRAC/IRS match- up, the source is not 
obvious, since both instruments require an extended source correction factor. A further consideration in our work is 
that we must apply a correction factor consistently across either all of the IRS or all of the IRAC data, since we only 
have overlapping information in a few regions, but we aim to determine the PAH fractions over the whole S'^MC map. 
There is no obvious way to decide which to correct, the IRS or IRAC calibration, so we choose the simplest route: 
we apply corrections to the IRS spectra based on the factors necessary to match IRAC 8 /im and MIPS 24 /im. We 
determine a correction factor that depends linearly on wavelength which matches the photometric to spectroscopic 
points in each region. This correction factor is very small in the region of the SL1/LL2 overlap, so it does not affect 
our stitching of the orders. 
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